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ABSTRACT 


Water is a key ingredient for the emergence of life as we know it. Yet, its destruction and reformation 
in space remains unprobed in warm gas. Here, we detect the hydroxyl radical (OH) emission from a 
planet-forming disk exposed to external far-ultraviolet (FUV) radiation with the James Webb Space Telescope. 
The observations are confronted with the results of quantum dynamical calculations. The highly excited OH 
infrared rotational lines are the tell-tale signs of H20 destruction by FUV. The OH infrared ro-vibrational 
lines are attributed to chemical excitation via the key reaction O-H;— OH-«H which seeds the formation of 
water in the gas-phase. We infer that the equivalent of the Earth ocean's worth of water is destroyed per 
month and replenished. These results show that under warm and irradiated conditions water is destroyed 
and efficiently reformed via gas-phase reactions. This process, assisted by diffusive transport, could reduce 
the HDO/H5O ratio in the warm regions of planet-forming disks. 


Water is a key ingredient in the emergence of life and is, therefore, a key aspect in the assessment of 
the habitability of (exo)planets. Yet, the trail of water to planets remains unclear. In a disk's inner regions 
(S 10 au), where the terrestrial and sub-Neptune planets are expected to form, a significant fraction of the 
water inherited from cold clouds is destroyed and reformed!—. In fact, theoretical models predict that 
under warm (T > 300 K) and FUV (6 « hv « 13.6 eV) irradiated conditions present in disk's atmospheres, 
water is destroyed and reformed by the chemical cycle'* O = OH = H20. Notably, this cycle is one of 
the processes that could lower the deuterium-to-hydrogen enrichment of water inherited from cold clouds 
and explain the intermediate D/H ratio found in Earth's oceans‘. 

In the past, understanding the detailed balance between the formation and destruction processes in a 
protoplanetary disk setting has been challenging because of the limited observational information. Warm 
water has been detected in planet-forming disks with Spitzer-IRS>-’ and ground-based instruments? but 
its ongoing chemistry remains unprobed. Here, we unveil the water cycle under warm and irradiated 
conditions thanks to the combination of JWST observations of OH and quantum dynamical calculations. 

The observational data stem from the JWST Early Release Science (ERS) program “PDRs4All””, 
which performed spectroscopic observations of the Orion Bar, an interstellar cloud exposed to the intense 
FUV radiation from the massive stars of the Trapezium cluster! (FUV flux at the ionization front equal to 
2—7 x 10^ times that of the local interstellar medium!! (ISM)). The Solar System presumably formed 
in a similar environment in which external UV radiation influenced the temperature, mass budget, and 
chemical composition of the solar nebula!?. Fig. 1 shows the d203-506 disk, located in the observed area, 


as a dark lane against a bright background of the nebulae (see Methods for its detailed characteristics). A 
photoevaporative wind, seen as a more diffuse emission of H5, is launched from the upper layers of the 
disk'?. This is the result of the intense FUV radiation from the massive stars in the proximity of d203-506 
that heats the gaseous disk. In addition, the young star at the center launches a collimated high-velocity 
jet, seen in [Fell] line emission which carves out a cavity in the wind. The North-West part of the cavity is 
particularly bright in excited molecular lines. This enhanced emission is likely due to a local increase in 
density induced by jet bow-shocks!^, and the direct irradiation from the Trapezium stars or the central 
star itself. Here, we report the detection of infrared rotational and ro-vibrational emission of OH by both 
MIRI-MRS and NIRSpec, respectively. The OH emission is detected where there is bright H? emission, 
namely all over the photoevaporative wind, revealing active gas-phase oxygen chemistry fueled by warm 
Ho (see also Extended Data Fig. 5). In this study, we focus on the bright spot to obtain high signal-to-noise 
Spectra. 


Evidence for H;O photodissociation. The MIRI-MRS spectrum shown in Fig. 2 reveals a series of 
highly excited rotational OH lines in the ground vibrational state, corresponding to a change in rotational 
quantum number of N — N-1. In total, lines from 23 rotational levels from N = 44 down to N = 18 are 
detected in the 9 — 13 um region, which probes upper-level energies as high as 45,000 K (see Extended 
Data Fig. 2). The excitation diagram of OH reveals an extremely high excitation temperature of 10,000 
K. This is ten times larger than the gas temperature inferred from the Hp rotational lines (see Extended 
Data Fig. 3 and 4). The detection of extremely rotationally excited OH is the smoking gun of water 
photodissociation. As already acknowledged from previous Spitzer observations?-!7 based on seminal 
molecular physics studies!* !?, the only process that can excite these lines is HyO photodissociation, via 
its B electronic state by short wavelength FUV photons forming OH in high-N states. The levels that 
are directly populated by water photodissociation (N ~ 35 — 45) were not detected before due to the 
limited spectral resolution and sensitivity of Spitzer-IRS. Furthermore, the higher spectral resolution of 
MIRI-MRS allows us to separate the symmetric from the anti-symmetric component of each rotational 
line down to 9.4 um (see Fig. 2 and Extended Data Fig. 2). Compared to the less excited lines observed 
with Spitzer!” , we find a much higher ratio (>10) between the symmetric and anti-symmetric lines. This is 
consistent with recent quantum calculations using a full-dimensional wave packet method which attributes 
this propensity to the Ë — X conical intersection pathway. 

In order to analyze the MIRI-MRS spectrum, synthetic JWST spectra were computed using the 
GROSBETA code?!, which takes into account the production of OH in excited states (see Method). As an 
input, we use the state distribution of nascent OH produced by H5O photodissociation via its first two 
electronic states as computed by ref.?:??, As shown in Fig. 2, the model agrees well with the observed 
spectrum. In particular, we recover the steep increase in line intensity from A = 9.15 um onward, and the 
relatively constant line intensities beyond 9.5 um. This behavior is intrinsically related to the rotational 
distribution of nascent OH (see Extended Data Fig. 6). Once an OH molecule is formed in a rotationally 
excited state, it decays via the N — N — 1 radiative transitions, a process called "radiative cascade "?!. The 
agreement between the model and the observations is not only an astronomical confirmation of a basic 
molecular process but, as shown in the next section, also provides access to the destruction rate of water 
and its local abundance. 


Evidence for chemical ("formation")-pumping by O+H,. The v = 1 — 0 ro-vibrational lines of OH 
are detected at shorter near-infrared (near-IR) wavelengths by NIRSpec, with upper rotational quantum 
number from N = 1 up to 10 (see Fig. 3 and Extended Data Fig. 2). These transitions have been detected 
in very dense environments like the innermost regions of protoplanetary disks? (ng > 108 cm~3), but 
not in a lower density environment such as the outer layers of the d203-506 system. The excitation of 
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the rotational levels within the v — 1 state is well described by a single excitation temperature of about 
1,000 K (Extended Data Fig. 4), close to the gas temperature inferred from H» emission (Extended Data 
Fig. 3) and ten times smaller than the excitation temperature of the mid-IR lines, implying a different 
origin. 

We attribute the near-IR OH lines to formation-pumping via the reaction O+Hz2, which is known to 
produce OH in vibrationally excited states”, but has hitherto never been observed in space. We amended 
the GROSBETA model to include the excitation of OH via the O--H» reaction using the state distribution 
of nascent OH extracted from the recent quantum calculations of ref.?^ (see Method and Extended Data 
Fig. 6). Since the state distribution of nascent OH depends on both the gas temperature and the state 
distribution of H5(v, J), we used the observed lines of Hz to compute, self-consistently, the distribution of 
nascent OH. The synthetic GROSBETA model shows good agreement with the NIRSpec spectrum (see Fig. 
3). The OH v = 2 — 1 lines are also predicted to be an order of magnitude weaker, in line with the upper 
limit set by the observations; quantum calculations predict that 2096 of the OH is produced in the v — 1 
state, with less than 1% in higher excited vibrational states (v > 2). The emission process is relatively 
simple, as an OH product formed in a v — 1 state rapidly cascades down via the ro-vibrational transitions. 
Other processes are known to produce OH in vibrational states (UV pumping, H20 photodissociation via 
its À electronic state, IR pumping) but we find that only inelastic collisions could significantly contribute 
to the observed emission if OH is very abundant (x(OH) > 1075, see detailed discussion in Method). The 
good agreement between observations and model points to formation-pumping of OH, thereby supporting 
earlier proposal based on Spitzer-IRS and Herschel observations! ^??. 


Warm gas-phase oxygen chemistry in action Because the newly formed OH products immediately 
decay via the observed transitions, the fluxes of these lines are directly proportional to the formation rate 
of OH via either route (see Table 1, and Method for further details). From the mid-IR lines, we derive an 
amount of water photodissociated per unit of time of  — 1.1 x 10? cm~? s~!. Considering the size of 
the bright spot (1.53 x 10^? sr, at a distance of 414 pc?), this corresponds to, roughly, the equivalent of 
the Earth ocean's worth of water being photodestroyed per month!. Similarly, from the near-IR lines, we 
infer an amount of OH formed via O-H» per unit of time of R = 1.2 x 10!! cm? s7!. 

These formation rates shed new light on the water cycle under warm and irradiated conditions. The 
mid-IR lines of OH reveal that water is actively photodestroyed in the d203-506 system. Water lines 
are not detected with MIRI-MRS, as predicted by non-LTE models (see Extended Data Fig. 8), likely 
due to the low density of the gas which is much lower than the critical densities of the H5O transitions 
observable with JWST. Yet, from the amount of H2O photodissociated per second, we infer a substantial 
amount of unseen H20 in the warm gas with N(H20) ~ 1—5 x 10? cm? (see Method). This quantity is 
still two orders of magnitude too low for water UV-shielding in this environment’. As estimated from 
the local FUV radiation field, the destruction timescale of water is about a day, much shorter than the 
dynamical timescale of the outflowing gas, which is longer than half a century (see Method). Therefore, 
the illuminated water cannot originate from inherited water formed in cold conditions. In fact, our analysis 
of the near-IR lines of OH demonstrates that the formation rate of OH via O+H)? is an order of magnitude 
larger than the destruction rate of H20. Water is thus efficiently replenished by the gas-phase formation 
route initiated by the reaction O-H» — OH+H. Our estimates also indicate that only a small fraction of 
OH, about P/R ~ 10%, is effectively converted into H20. Notably, OH is also photodissociated, reducing 
the probability of OH being converted into H5O. Interestingly, OH photodissociation has already been 
unveiled by the Hubble Space Telescope observations of [OI] emission at A = 6300 À?5, which is likely 
excited by OH photodissociation??. Assuming that [OI] emission is produced via OH+hv — O+H with a 


'The oceans of the Earth are composed of about 5 x 1076 molecules of water. 
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50% probability to form oxygen in the O(! D) state??, we infer an amount of OH photodissociated per unit 
of time of 7x10? cm"? s~!; a value close to the formation rate R via O--H» given above. 

What are the implications for the trail of water to terrestrial planets and notably to Earth? 
The enhanced D/H ratio in standard mean ocean water (SMOW)?! of 1.5 x 1074 relative to the bulk 
interstellar elemental D/H ratio (~ 2 x 10~>)** indicates that a fraction of terrestrial water formed under 
cold conditions, likely at the surface of interstellar grains??. This anomaly would indeed reflect the effects 
of chemistry at low temperatures where the small zero-point energy difference between D- and H-bearing 
species can create large deuterium fractionations*+ 35. However, the observed D/H ratio in protostars*°-** 
— tracing the inherited water content — is higher than the D/H ratio in SMOW (0.3 — 1 x 1073 versus 
1.5 x 107^, respectively). Hence, chemical processing could have occurred in warm gas, reducing the 
deuterium fractionation because the formation rates of HDO and H5O are similar at high temperatures. 
In the d203-506 system, H5O is efficiently destroyed and reformed via the warm gas-phase formation 
route O — OH =€ H20. This can significantly decrease the water D/H ratio inherited from the cold phase 
to interstellar elemental values (~ 2 x 107?)??. The resulting HDO/H20 ratio could be even lower if 
selective photodissociation of dihydrogen leads to subsolar values of HD/H»??. It is however unclear if 
water unveiled in the d203-506 system is incorporated into comets and asteroids since the observed OH 
might rather trace unbounded gas heated by the external FUV radiation. The latter is further supported by 
models of photoevaporative winds with similar conditions as d203-506, which predicts photodissociated 
OH (and consequently H20) to be primarily present in the unbounded gas ^. 

Nevertheless, the d203-506 system constitutes a unique interstellar laboratory as the same chemical 
processes will occur in the atmosphere of planet-forming disks even if the disk is solely irradiated by 
the host star. The FUV fields generated by accreting stars at the disk surface range from Go —10? to 
10° between 1 and 10 au, depending on the dust properties and shielding by the gas^?"; a range that 
is comparable to that of the gas probed by our observations. The local density inferred from OH mid- 
IR lines (see Method) is somewhat lower compared with gas densities expected in disk upper layers 
(0.2 — 1.2 x 107 cm? versus 105 — 101? cm?) and the temperature is higher than that found for water in 
planet-forming disks (1000 K versus 400-600 K)’. Still, the warm gas-phase formation route is already 
efficient at T > 400 K and the gas density only affects the relative OH/H5O ratio in the irradiated layers?! . 

We expect mid- and near-IR emission of OH to be detected by JWST in a wide range of astrophysical 
environments (e.g., shocks, photodissociation regions at the edge of molecular clouds, comets). Our new 
analysis method of OH lines will provide vital constraints on the physics and chemistry of the gas in such 
environments. This work also constitutes a benchmark illustration of the potential of excited IR lines 
from small molecular species, which can be transposed to other species like CO, CH+42 or CHj. In this 
context, new quantum dynamical calculations and laboratory experiments are essential to exploit the full 
potential of JWST. 
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Jet 


Figure 1. JWST NIRCam composite image of the Orion Bar, located in the Orion molecular cloud. Red is the 3.35um emission (F335M NIRCam filter), blue 
is the emission of Paw (F187N filter subtracted by F182M filter) and green is the emission of the Hz 0-0 S(9) line at 4.70 um (F470N filter subtracted by 
F480M filter). The inset shows a zoom-in of the d203-506 planet-forming disk where OH lines are detected. Red is the emission of the Hz 1-0 S(1) line at 
2.12 um (F212N filter), blue is the [Fel] line emission at 1.64 um (F164N filter), and green is the emission in the F140M broad-band near-IR filter that traces 
scattered light around 1.4 um (F140M broad band filter). The white contours represent the emission of the OH rotational line at 9.79 um detected with 
MIRI-MRS (levels are 1.1, 2.5 x 10? erg cm~? s^! sr-!). The bright spot in the northwestern part of the d203-506 system coincides with the region of 
interaction between a jet and the photoevaporative wind. The brightest OH emission originates from here. The spectra shown in Fig. 2 and 3 are an average 
over the region delineated by the green circle in order to have the best S/N. Image credits: NASA, ESA, CSA, PDRs4All ERS Team, pdrs4all.org. 
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Figure 2. Evidence for H20 photodissociation from OH rotational lines detected by MIRI-MRS. (Top panel) Observed MIRI-MRS spectrum where the dust 
continuum and bright lines, other than that from OH, have been subtracted (see Method). Rotationally excited OH lines are detected up to levels N — 44. 
(Bottom panel) Synthetic spectra from GROSBETA. The spectrum in red assumes that water photodissociation produces OH in symmetric A-doubling states, 
whereas the gray spectrum assumes equal distribution among the sub-levels. 
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Figure 3. Evidence for formation-pumping via O+H2 from OH ro-vibrational lines detected by NIRSpec. (Top panel) Observed NIRSpec spectrum where the 
dust continuum and the bright lines, other than those from OH, have been subtracted (see Method). Each ro-vibrational transition is split by the spin-orbit 
coupling, however the A-doubling is unresolved. (Bottom panel) Synthetic spectra from GROSBETA including the excitation of OH via formation-pumping. 
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Table 1. 


Parameters derived from the analysis of OH and H; lines 


Diagnostics Quantity Measured value 

H; lines N(H3) (9.9+0.9) x 10? cm? 
T° 950 +17 K 

OH mid-IR lines | ® (rate HO Av)? | (1.1 £0.2) x 109 cm? s7! 
N(H20) 1 x10!5 - 5 x 10!5 cm~? 
x(H20)° 5x 1076 -2.5 x 1075 

OH near-IR lines | R (rate O-H5)^4 (1.20.3) x 10! cm? s7! 
ny (1.43 0.4) x 107 cm? 


inferred from Hp lines. 


s land ny = 1.310° cm 3. 


Note a: the uncertainties of these values do not take into account calibration effects. They could be 
increased by a factor of 2-3 if the uncertainty associated with calibration effects is equal to 20%. 

Note b: number of OH molecules formed per unit area and time via either formation route (in cm 
Note c: the abundance with respect to the total number of hydrogen atoms is calculated using N(H2) 


"uc 


Note d: the inferred value would be reduced if inelastic collisional excitation were to play a significant 
role. The lower limit expected for these values are R = 1.1x10!? cm? 


3 
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Methods 


Observations and data reduction Observations with MIRI-MRS in the integral field unit (FU) mode 
were obtained as part of the Early Research Science (ERS) program PDRs4All: Radiative feedback from 
massive stars (ID1288, PIs: Berné, Habart, Peeters)?. The observations cover a 9 x 1 mosaic centered 
on 042900 = 05'357in20.47495, dy2009 = -05?25' 10.45". A 4-point dither optimized for extended sources 
was applied and the FASTR1 readout pattern adapted for bright sources was used. The integration time 
was 521.7s using 47 groups/integration and 4 integrations. The data span a wavelength range from 4.90 
to 27.9 um, have a spectral resolution R ~ 1700-3700, and a spatial resolution of 0.2"-0.8". The latter 
corresponds to a very small spatial scale of about 100 AU at the distance of the Orion Bar (414 pc?9). 
The MIRI-MRS data were reduced using version 1.11.1 of the JWST pipeline”, and JWST Calibration 
Reference Data System? (CRDS) context 1097. The stage 2 residual fringe correction was applied in 
addition to the standard fringe correction step. A master background subtraction was applied in stage 3 
of the reduction. The 12 cubes (4 channels of 3 sub-bands each), all pointing positions combined, were 
stitched into a single cube (Canin et al., in prep.) (see ref.*” for data reduction details). The mosaic was 
positioned to overlap with the PDRsA4AII NIRSpec IFU observations of the Orion Bar which is also a 9 
x | mosaic. The NRSRAPID readout mode appropriate for bright sources, and a 4-point dither pattern 
were also used for NIRSpec. The on-source integration time was 257.7 per exposure with five groups per 
integration with one integration per exposure. The data were reduced using the JWST Science Calibration 
Pipeline (version 1.10.2.dev26--g8f690fdc) and the context jwst_1084.pmap of the Calibration References 
Data System (CRDS) (see ref.!! on data reduction process). The MIRI-MRS and NIRSpec spectrum, 
in units of MJy sr !, were averaged in a circular aperture centered on Œy2000 = 5h35min50 31455 and 
6120007-5^25' 05.528" with a radius of 0.15". An OFF position was averaged in a circular aperture centered 
on 090p = 5°35™"20.25398 and 5j2999=-5°25’ 05.498" with a radius of 0.3". 


Post-processing of the JWST spectra. In order to better visualize OH lines in Figs. 2 and 3, the 
continuum was subtracted considering several points in the spectral regions free of line emission and 
lines other than that of OH were removed as shown in Extended Data Fig. 1. In the selected MIRI-MRS 
spectral range (9-10.4um), several prominent lines are present, such as HI lines 8-13 at 9.261 and 9-20 at 
9.392 um, H? 0-0 S(3) line at 9.665 um. The lines from the background of the planet-forming disk were 
subtracted by fitting them in the OFF spectrum. This allows subtracting the emission of the recombination 
lines of HI, 8-13 at 9.261 and 9-20 at 9.392 um. We also subtracted the H? 0-0 S(3) line by fitting it 
directly in the ON spectrum (see Extended Data Fig. 1). The same processing was applied to the NIRSpec 
spectrum; a continuum was subtracted and HI lines, 5-11 at 2.873 um and 5-10 3.039 um, Ol lines at 
2.893 and 3.099 um, H» lines 2-1 O(3) at 2.974 um, 1-0 O(4) at 3.004 um, 2-1 O(4) at 3.189 um and two 
unidentified lines at 3.164 and 3.224 um were removed by fitting them directly in the ON spectrum. 


Planet-forming disk d203-506 characteristics. The planet-forming disk is an almost edge-on disk seen 
in silhouette against the bright background. The measured radius is Rout = 98 + 2 au and the total mass is 
estimated to be about 10 times the mass of Jupiter*®. Kinematic studies with ALMA point to a stellar 
mass of the host star to be below 0.3 Mo. This environment is the archetype of a young planet-forming 
disk in the intermediate phase between being embedded and being unshielded in the surrounding HI gas 
irradiated by an external FUV field. d203-506 shows no signs of the presence of an ionization front^* ??:45, 
indicating that the radiation field reaching the disk is completely dominated by FUV photons (E « 13.6 
eV). The Far-UV (FUV) radiation field incident on the ionization front (IF) of the bar is Go = 2 — 7 x 10* 


?https://jwst-pipeline.readthedocs.io/en/latest/ 
5https://jwst-crds.stsci.edu/ 
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in Habing units! ^^? (Go = 1 corresponding to a flux integrated between 91.2 and 240 nm of 1.6 x 10? 
erg cm? s7!) as derived from UV-pumped IR-fluorescent lines! ^9, d203-506 can be illuminated by 
the Trapezium cluster (©! Ori C) but also by 8? Ori A?!. The estimation of the FUV field at the surface 
of the disk is uncertain but it is expected to be similar than at the IF as determined with measures by 
geometrical considerations (distance between the disk and both stars) and UV-pumped IR-fluorescent 
lines (OI and H5)?. Moreover, Hz ro-vibrational lines highlight a bright molecular emission originating 
from a photevaporative wind surrounding the edge-on disk. As seen in Fig. 1, there is a bright emission 
spot visible in the Hz emission in the northwestern part of the d203-506 system. This spot appears to 
coincide with the region of interaction between a jet and the wind. It is only visible on the side facing the 
Trapezium. The molecular lines (such as OH or pure rotational H» lines) between the bright spot and the 
more diffuse wind are similar in relative intensity, suggesting that the molecular excitation is primarily 
powered by the external FUV. However, shocks could impact the bright spot's chemistry because the 
collimated jet could compress the gas, increasing the density locally, explaining the brightness of the spot. 

In this bright spot, we detect a very active warm water chemistry. We can estimate two timescales to 
understand its origin. The chemical timescale for water chemistry is determined using the photodissociation 
rate of water kp: Tphoto = | /kg. This estimation leads to a chemical timescale of about a day. In 
comparison, the dynamical timescale of the outflowing gas Tayn is much longer. It can be estimated as 
Tayn = d /Vaay, Where d is the distance between the inner disk and the bright spot and vagy is the advection 
velocity. In the case of d203-506 d ~ 100 au and Tayn Z 50 years. 


Excitation diagram of H; and OH. The integrated line intensities of OH and H? measured with NIRSpec 
and MIRI-MRS are reported in Extended Data Tables | and 2. They were derived by fitting each line in the 
spectra by a Gaussian function plus linear function to account for the continuum. We note that the H2 line 
intensities toward d203-506 are also reported in ref.?? using a slightly different aperture. The excitation 


Az 
hvoAuigu 


Eup of each transition, where / is the integrated line intensity, and A,,; and Vo are the Einstein-A coefficient 
and the frequency of the line, respectively. An excitation temperature Tọ, and, for H» lines, a total column 
density N can then be inferred using the canonical formula 


diagrams of OH and H} are derived by plotting Y = In ( ) as a function of the upper energy level 


N 
Y=In (sa) — Ey/KpTex, (1) 


where Q(Tz) is the partition function. 

The excitation diagram of OH lines is presented in Extended Data Fig. 4 and analyzed in the main text. 
The excitation diagram of H» shown in Extended Data Fig. 3 was made using the PhotoDissociation 
Region Toolbox Python module? (pdrtpy)??-76. The excitation of the ro-vibrational lines up 
to energy levels of ~ 8,000 K can be described by a single excitation temperature of ~ 1,000 K. This 
corresponds to the average gas temperature of the warm H» layer where OH and H20 are efficiently 
formed and photodestroyed. We further infer from Eq. (1) a column density of N(H2) = 9.9 x 10? cm ?. 


The OH model. Our OH model includes the ro-vibrational, A-doubling, and fine structure levels of 
OH in its ground electronic state OH(X?IT) and in its first electronic state OH(A?X*) as provided by 
refs.??^:?^ and compiled in ref.?!. The synthetic JWST spectra of OH are calculated using the molecular 
excitation code GROSBETA”! and the results of quantum dynamical calculations!???, The excitation 
model is based on a single-zone approach following the formalism presented in ref.??. The population 
of OH levels is computed by considering the production of OH in various states, radiative pumping, and 


^https://github.com/mpound/pdrtpy 
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inelastic collisional (de-)excitation. Assuming chemical steady-state, the detailed balance equation for 
level i along the line of sight is given by: 


(2) 


iud ay 


Y PFjNj-NiY PF x (nom - NOH 


jfi jfi 
where N(OH) is the total column density of OH and N; [cm ?] are the column densities in level i. P;; are 
the radiative and inelastic collisional transition probabilities. The inelastic collisional (de-)excitation of 
OH with He and H, was determined using inelastic collisional rate coefficients from ref.?9?" and have 
been extrapolated to include collisional transitions between higher rotational levels of OH such as done 
in ref.?!. In this work, the column density of OH is fixed at N(OH) = 2 x 10? cm"? and the density at 
ng = 107 cm^?. We further adopt a local radiation field composed of a diluted black-body at 40,000 K to 
model the UV radiation field emitted by the Trapezium stars, an infrared radiation field corresponding to 
that detected by NIRSpec and MIRI-MRS, and a diluted black-body 50 K to account for the far-IR and 
mm emission of the dust. We stress that because inelastic collisional de-excitation and radiative pumping 
of the OH levels probed by JWST is orders of magnitude weaker than radiative decay, the exact values of 
N(OH), ny, and local radiation field, are of little importance. 

In Eq. (2), F is the formation rate, measured in cm ? s~!, that produces OH in specific states i with a 
probability denoted as f; and we assume that the probability of destroying OH from a state i is equal to the 
proportion of OH in this state N;/N (OH). We assume that the other formation and destruction processes do 
not impact the level population. For the sake of clarity, we compute synthetic spectra including either H20 
photodissociation or chemical formation-pumping. This approach is valid because these two processes 
excite different energy levels. Thus, we define two different formation rates: the formation rate via the 
photodissociation of water denoted as ® and the formation rate via O+H3 denoted as R. ® is also the 
column density of H20 photodissociated per unit of time?!: 


pz f konmoaz, (3) 
zZ 


where np,o is the local number density of H20 ([cm ?]) and ko is the photodissociation rate of H2O 
forming OH ([s~!]), which depends on the strength and shape of the local FUV radiation field^?. Similarly, 
R is also the column density of OH formed per unit of time via O+H). 


R= [ Y X kii(T)x;(H)nmnodz= [ imi nodz, (4) 
Zi j z 


where no and ng, are the number densities of atomic oxygen and molecular hydrogen, x;(H») is the level 
population of H», and k;_,; is the state-to-state rate coefficient of the reaction ([cm? s-!]). 

The distributions of nascent OH /; stem from the results of quantum dynamical calculations and are 
shown in Extended Data Fig. 6. For H2O photodissociation, we use the rotational and vibrational state 
distribution of OH produced by H20 photodissociation via both its A and B electronic states calculated by 
ref.?? and ref.!° and compiled by ref.?!. Since these pioneering calculations did not include the spin-orbit 
coupling nor the A-doubling, we include the results from ref.?? by assuming that photodissociation via the 
B state leads to OH in symmetric states, with equal distribution between the two spin-orbit manifold (see 
Extended Data Fig. 2 for the spectroscopy of OH). The resulting state distribution of OH is calculated by 
integrating the state-specific cross-section over a UV field represented by a black-body at Te sf = 40,000 K 
(see ref.?!). In the absence of systematic quantum dynamical calculations for short wavelength photons, 
we further neglect the impact of H20 photodissociation shortward of 114 nm on the excitation of OH. 
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For formation-pumping of OH, we extracted the distribution of nascent OH from the state-to-state rate 
coefficients computed by Veselinova et al. in ref.?^. Because the distribution of OH depends on both the 
kinetic temperature and the population of H», the distributions are computed as 


fi(OH) = ) kj oi(T)x;(H2)/ ki (T); (Ho), (5) 
J E, 


where x;(H») is the level population of H2 and &;.,; (T) are the state-to-state rate coefficients at temperature 
T. The values of x;(H2) are inferred from the excitation diagram of Hz, which probes the majority of 
Hp levels. For the Ho(v = 0,J = 0), H2(v = 0,J = 1) and H2(v = 0,J = 2) levels, which are not probed 
by JWST, we extrapolated their populations assuming local thermodynamical equilibrium. Interestingly, 
we find that the resulting distribution of OH differs very little from a distribution assuming a Boltzmann 
distribution of H» levels at T ~ 1,000 K. In fact, the formation rate of OH(v = 1) is driven by the 
rotationally excited levels of Hz (v=0, J ~ 4 — 9) that are much more populated than the vibrational levels 
and are reactive enough to produce vibrationally excited OH by reacting with O atoms. 


Conversion between line intensity and formation rate In the d203-506 system, we find that the highly 
excited rotational and vibrational OH levels are primarily populated by formation pumping (the reaction 
O+H; or the photodissociation of water) followed by radiative decay, with negligible impact of radiative 
pumping and inelastic collisions. Furthermore, the ro-vibrational and rotationally excited lines are found 
to be optically thin. In the optically thin regime, the integrated line intensity is proportional to the column 
density of OH in the upper energy state: 


lij = c Ni. (6) 


Assuming that excited OH levels are only populated by formation pumping and radiative decay and that 
the fraction of OH in rotationally and ro-vibrationally excited states constitutes a negligible fraction of the 
total population of OH, the detailed balance equation (2) simply writes: 


Mjc-Aj if izj 
Y MiNj+F x fi=0 with 4Mj—-YAs if i=j (7) 
jj£i k 


This system of linear equations on N; shows that N;, and therefore the intensity J;;, are simply proportional 
to F, the formation rate of OH. Following ref.?!, we can therefore define the intensity as 


hV; j 
An 


lij — l;jF, (8) 
where /; j is the dimensionless proportionality factor. It corresponds to the probability that a nascent OH 
eventually cascades via the radiative transition i — j. 

I; j depends only on f;, which in turn depends on the shape of the radiation field when dealing with 
excitation by H5O photodissociation, and the gas temperature and the population of H5 when dealing with 
chemical formation pumping. For water photodissociation, /; j 18 provided in Appendix D of ref.?! for 
various shapes of the UV radiation field. For the chemical pumping formation via O+H2 — OH-+H, the 
proportionality factor computed at T = 1,000 K summed over all the four component of the v = 1 — 0, 
N = 4 — 5 ro-vibrational lines is / = 3 x 107°. 
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Estimation of column density of warm water and local gas density Our knowledge of the evolution 
of matter in space, from diffuse clouds to planets, depends on our ability to assess the local physical and 
chemical properties. Thanks to the derived value of the OH formation rates via the photodissociation of 
water ® and via O+H2 R, one can derive astrophysical quantities such as the column density of warm 
water N(H5O) and the local gas density. Assuming an homogeneous medium Eq. (3) can be rewritten as 


o 


where the photodissociation rate kg depends on the FUV radiation field. Based on the OH mid-IR lines, 
we estimate ®=1.1 x 10!° cm? s^! (see main text). We can obtain an estimate of the photodissociation 
rate of H20 based on the local FUV field intensity. In the d203-506 system, the exposed reservoir of water 
likely lies at low visual extinction. Here, we assume that the FUV field at the H?/H; transition is extincted 
by as much as a factor 5 such that kg ~ 1 x 107? — 2 x 1076 s- in the bright spot, corresponding to a 
visual extinction of Ay < 1. This translates to N(H50)- 1 x 10/5 — 5 x 10P cm-?. 


Similarly, assuming a homogeneous medium, the local density can be estimated from the inferred 
value of R from Eq. (4) as: 


R 
k(T)N(H2)x(0) 


ny = (10) 


where k(T) is the rate coefficient of the reaction, N(H2) is the column density of H2, x(O) is the oxygen 
abundance, and ny is the total number density of hydrogen atoms. For d203-506, we found a column 
density of warm excited Hy of N(H2) = 9.9 x 101? cm~? and T ~ 1,000 K from the rotational and 
ro-vibrational lines (see Extended Data Fig. 3). From the estimation of temperature and excitation of H2, 


we can derive the coefficient rate k = 2.8 x 1071? cm? s~!. Using the estimate of R = 1.2x10!! cm~! s7! 
g 


from the OH near-IR lines, we derive ng = 1.4x 107 cm 3. 

Ro-vibrational excitation of OH In this work, we provide evidence that the near-IR OH lines are excited 
by chemical pumping via O+H, — OH-H. In this section, we review alternative processes that can lead to 
vibrationally excited OH in the ground electronic state. 

Water photodissociation via the À electronic state by long wavelength photons produces rotationally 
cold, but vibrationally hot OH. Our OH model does include this process using the OH distribution 
computed as described in ref.!? while using the improved H20 A state potential energy surface of ref.??. 
The resulting OH emission spectrum when considering only water photodissociation via both the À and 
B states is shown in Extended Data Fig. 7. The overall shape of the modeled near-IR spectrum is very 
different from that observed. The ro-vibrational line within the v=1-0 band exhibits a much lower rotational 
temperature, and the v=2-1 band is predicted to be strong. Adopting the value of d that fits the mid-IR 
lines of OH, we also find that the ro-vibrational lines typically a factor 10-20 weaker than the observed 
near-IR lines. 

UV radiative pumping of OH via the A-X band can also produce vibrationally excited OH in the ground 
electronic state. Our excitation model does take into account this process since the OH levels extend to the 
non-dissociative OH(A?X-^) electronic state. Discarding OH chemical pumping in the excitation model, 
we find that the OH near-IR spectrum is about two to three order of magnitude weaker than observed, 
and exhibits a too-cold rotational temperature within the v = 1 state. We note that for UV pumping, the 
intensity of the lines depends on the assumed column density of OH, which is not directly probed by 
our observations. However, the adopted column density of N(OH) = 2 x 10? cm? is a good first-order 
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estimate, since it corresponds to an OH abundance of x(OH) ~ 1075, a value that is in line with the results 
from physicochemical models*!. 

Little is known about the inelastic collisional cross sections connecting the v = 0 and the v = | states 
of OH. Owing to the repulsive barrier between OH and H»5, and OH and He, collisions with H» and He are 
negligible compared to those with atomic hydrogen, which is likely abundant in the OH emitting layer. 
Collision with the latter is thus the only process that can contribute to the excitation of OH(v = 1). In 
fact, ref. computed the inelastic de-excitation collisional rate coefficients for OH(v=1)+H — OH(v=0) + 
H below a kinetic temperature of T = 300 K and assuming a thermal distribution of OH levels. In the 
d203-506 system, the temperature is much higher, about T ~ 1,000 K, and the rotational levels of OH 
even in the v = 0 are predicted to be subthermally populated due to the low gas density (ny ~ 10’ cm~°). 
Extrapolating the result of ref. to the conditions in the OH emitting gas remains uncertain. Still, 
for a temperature of T ~ 1,000 K, we estimate a rate of the order of kco ~ 5 x 107cm? s~!. The 
ratio between the OH(v=1) formed by chemical formation-pumping and by inelastic collisions is about 
x(OH)keon /x(O)kchem,v=1, Where kcnem,v—1 is the formation rate O--H» leading to OH(v = 1). We conclude 
that chemical pumping dominates over inelastic collisional excitation for an abundance ratio of about 
x(O)/x(OH) = 10. The calculation of inelastic collisional rate coefficients is warranted to evaluate the 
contribution of inelastic collisions in the emission of OH. If inelastic collisions were to contribute to the 
near-IR emission of OH, the inferred chemical formation rate R would be lower. However, we stress that 
the formation rate of O+H2 — OH+H is close to the destruction rate of OH through photodissociation, as 
inferred from [OI] emission. This provides further evidence that the bulk part of OH emission is produced 
via chemical formation and not by inelastic collisional excitation. 


Quantum dynamical calculations for O+H,. This section summarizes the ab initio quantum deter- 
mination of the state-to-state rate coefficients of the reaction O (P) + Ho (v, j) — OH (v/, j^) +H used in 
our OH excitation and chemical models, see also ref.?^. Time-independent quantum mechanical (QM) 
scattering calculations were carried out applying the coupled-channel hyperspherical method implemented 
in the ABC code?! yielding the usual S matrix, from which the state-to-state cross sections and rate 
coefficients can be computed. The set of potential energy surfaces (PESs) calculated by ref.°* was used, 
which includes the PESs of the two electronic states (one of symmetry ?A" and one of symmetry ?A' ) that 
correlate adiabatically with the electronic ground state of reactants, O^ P)-H», and products, OH(7T1)+H. 
Scattering calculations were carried out at 60 energies between 0.3-2.5 eV, including all partial waves (J) 
needed to reach convergence (Jmax = 62) and all the values of the helicity (Q) up to 26. QM scattering 
calculations provided the cross section as a function of the rotational quantum number assuming that 
OH is a closed-shell molecule. To calculate the nascent rotational population of OH in terms of the total 
angular momentum excluding spin, N, we plotted the calculated cross sections calculated as a function 
of the rotational energy, and assigned to the N state the cross section corresponding to the actual energy 
of that state. Calculations do not distinguish between the two spin-orbit manifolds, and the procedure 
of ref. was used to distribute the population of the product's rotational states into the two A-doublet 
components. The propagation was carried out in 300 log-derivative steps, with a maximum value of the 
hyperradius of 20 ao. The basis included all the diatomic energy levels up to 3.25 eV. 

Since reactions can occur on both ?A' and 3A” electronic states, the total rate constant of the reaction 
is obtained as a weighted average of the respective rates ky, and kav. The weights are obtained by 
establishing the correlation between the electronic states A' and A" and the spin-orbit levels of oxygen, 
and are thus related to their relative population. It is found that the 13A” state correlates with three of 
the five components of OCP), while the 1?A' state correlates with the other two components of OP) 
and one component of OC Pj). Under these considerations and assuming a thermal distribution of the 
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spin-orbit levels of the oxygen atom, the overall thermal rate coefficients are calculated as ref.?^: 


7 3kan(T) + |2 +exp ( — AR) ky (T) 


T)= 11 
) 5+3-exp(— 4&5) + exp ( — ^52) Tu 


where k(T) referrers to the desired (v, j  v', j’) state-to-state rate coefficient at a given temperature, ky (T) 
and kar (T) are the state-to-state rate coefficients on the A’ and A" PESs, respectively. AE; 227.708 K and 
AE=326.569 K are the energies of the oxygen atom spin-orbit levels ?P; and ?P», respectively, over the 
ground state °P). We should note that the ?P» contribution and two contributions of the ?P; do not appear 
in the numerator; this reflects the fact that these spin-orbit states are not reactive and will not contribute to 
the total rate constants of OH formation. 
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Extended Data Figure 1. (a) (Top panel) Spectrum observed with MIRI-MRS. The red line is the estimated continuum. (Middle panel) Continuum 
subtracted spectrum. The red Gaussians are the fits to lines other than OH. The blue Gaussians are the fits to lines from the OFF position which contaminate 
OH lines. (Bottom panel) Processed spectrum with the dust continuum and the bright lines other than OH subtracted. (b) (Top panel) Spectrum observed with 
NIRSpec. The red line is the estimated continuum. (Middle panel) Continuum subtracted spectrum. The red Gaussians are fits to lines other than OH. (Bottom 
panel) Processed spectrum with the dust continuum and the bright lines other than OH subtracted. 
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Extended Data Figure 2. Schematic view of OH spectroscopy relevant to /WST. (a) OH rotational and ro-vibrational energy levels. The red arrows are the 
detected pure rotational transitions with MIRI-MRS. The blue arrows are some of the detected ro-vibrational transitions with NIRSpec. (b) Zoom on the yellow 
box to reveal the splitting of a rotational level due to the spin-orbit coupling and the A-doubling. The two spin-orbit states are labeled by the Q quantum number 
and the A-doubling states are labeled by their € = e /f spectroscopic parity. The green and pink arrows are the transitions detected in the observations forming 
a quadruplet. The green arrows are the transitions arising from symmetric states and the pink arrows are the transitions arising from anti-symmetric states. 
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Extended Data Figure 3. Excitation diagram of Hp at the bright spot. The temperature fit was made on the first five pure rotational lines. The measured line 


fluxes are reported in Extended Data Table 2. 
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Extended Data Figure 4. Excitation diagram of ro-vibrational and pure rotational lines of OH. The blue crosses are the ro-vibrational transitions detected 
with NIRSpec (v=1-0) and the red crosses are the pure rotational transitions detected with MIRI (v=0-0). The measured line fluxes are reported in Extended 


Data Table 1. 
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Extended Data Figure 5. Spatial distribution of the emission of key lines. (Top Left) d203-506 image from the NIRCam F212N filter (in MJy sr—!), (Top 
Right) integrated intensities image from NIRSpec OH line at 2.934 um (in erg cm? s^! sr~!), (Bottom Left) integrated intensities image from MIRI-MRS 
OH line at 9.791 um (in erg cm~? s^! sr~!), (Bottom Right) integrated intensities image from Hubble Space Telescope (HST) OI at 0.63 um (in counts s^! )?$. 
Contours of NIRSpec H3 1-0 S(1) line are shown in white (at 7 x 1074, 2 x 10? erg cm ? s^! sr P). 
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Extended Data Figure 6. State distributions of nascent OH predicted by quantum dynamical calculations summed over the spin-orbit and A—doubling 
states. (Left) Distribution of nascent OH following its formation via O--H» at a temperature of T = 1,000 K and for a H» population as inferred for d203-506. 
(Right) Distribution of nascent OH following H20 photodissociation by an FUV field representative of the Orion Bar. Photodissociation by short wavelength 
photons A < 144 nm via B electronic state of water produces rotationally hot OH with rotational quantum numbers of N ~ 35 — 45 and photodissociation by 
longer wavelength photons via the A electronic state leads to rotationally cold but vibrationally hot OH. 
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Extended Data Figure 7. OH near-IR synthetic GROSBETA models for (Top) UV radiative pumping and (Bottom) H20 photodissociation via its A 
electronic state. None of these processes can account for the shape and strength of the observed OH ro-vibrational spectrum. 
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Extended Data Figure 8. MIRI-MRS observations versus a synthetic RADEX?? spectrum of the brightest H5O lines adopting the maximum amount of 
unseen water of N(H)0) = 5 x 10? cm? as inferred from mid-IR lines of OH, a temperature of T=1000K as inferred from H3 lines and a density of 

ng = 107 cm? as inferred for near-IR lines of OH. Inelastic collisional rate coefficients are from ref.* available on the LAMDA database®. We also assumed 
an IR background inferred from NIRSpec and MIRI-MRS observations. When calculating the line intensities, we assumed that the IR continuum background 
interacting with the gas is not along the same line of sight as the observations are taken as described in ref.?! . The latter assumption provides a strict upper limit 
on the line strength. Owing to the low gas density, H20 is subthermally excited, leading to undetectable lines. In the 5-7 um region, the MIRI-MRS spectrum 
is affected by residual fringe and possible contamination by CHi 52. A LTE model at T = 700K of CH} is overlayed in grey. In the 18-24um, the continuum 
increases greatly so the noise increases accordingly making the detection of H20 lines in this region impossible. The same model of H2O multiplied by a factor 
of 10 is overlayed in purple. 
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MIRI-MRS 


À (um) | N | Eup/kp (K) | 1 (X102? erg cm? s~! sr) 
9.181 | 44 | 43,343 0.91 € 0.12 
9.219 | 43 | 41776 1.02 + 0.07 
9.266 |42 | 40216 1.50 + 0.06 
9.325 | 41 38,663 1.38 + 0.05 
9.394 | 40| 37,120 1.64 £0.16 
9.475 | 39 | 35,588 1.30 + 0.09 
9.568 | 38 | 34,070 1.89 + 0.08 
9.791 | 36] 31,079 2.97 021 
9.903 |35 | 29,609 3.24 + 0.23 
10.069 | 34 | 28,160 3.89 + 0.30 
10.232 | 33 | 26,731 3.86 + 0.12 
10.411 | 32 | 25,325 3.48 + 0.18 
10.608 | 31 23,943 2.58 + 0.18 
10.825 | 30 | 22,586 4.44 + 027 
11.610 | 27 18,671 3.38 + 1.35 
11.615 | 27 18,701 2.24 + 0.87 
11.929 | 26 17,447 3.30 + 0.26 
12.657 | 24 15,053 1.12 £0.10 
12.665 | 24 15,085 0.92 + 0.05 
13.081 | 23 13,916 1.27 £0.09 
13.089 | 23 13,949 2.20 + 0.18 
13.550 | 22 12,816 1.36 + 0.09 
13.560 | 22 12,850 1.24 + 0.12 
14.070 | 21 11,754 0.94 + 0.07 
14.082 | 21 11,789 0.84 + 0.04 
14.648 | 20 10,783 1.29 + 0.04 
14.664 | 20 10,767 1.27 £0.07 
15.295 | 19 9,749 3.34 + 0.44 
15.314 | 19 9,786 1.71 £0.08 
16.021 | 18 8,809 1.26 + 0.03 
16.045 | 18 8,847 1.16 £0.04 

NIRSpec 
A (um) | Q | N | Eup/ke (K) | I (x10? erg cm? s~! sr) 
2.718 | 1/2 | 3 5414 2.06 + 0.65 
2.851 | 1/2} 1 5314 2.89 + 1.02 
2.870 | 3/2} 1 5134 6.77 + 1.24 
2.886 | 1/2 | 2 5401 4.81 + 0.97 
2.901 | 3/2 | 2 5250 7.92 + 0.86 
2.922 | 1/2 | 3 5541 6.35 + 0.96 
2934 | 3/2 | 3 5415 10.00 + 0.41 
2.960 | 1/2 | 4 5735 4.99 + 0.81 
2.970 | 3/2 | 4 5627 7.99 + 0.61 
3.000 | 1/2} 5 5982 5.12 + 0.55 
3.008 | 3/2} 5 5888 10.18 + 0.98 
3.048 | 3/2 | 6 6199 7.78 + 0.31 
3.085 | 1/2] 7 6632 4.88 + 0.78 
3.091 | 3/2] 7 6558 7.28 + 1.17 
3.131 | 1/2] 8 7033 3.50 + 0.51 
3136 | 3/2] 8 6966 6.69 + 1.21 
3.179 |12|9 7483 4.18 + 0.64 
3183 | 3/2 9 7422 4.54 + 0.83 
3229 | 1/2 | 10 7483 2.95 + 0.93 


Extended Data Table 1. (Top) Intensities of the OH lines detected with MIRI-MRS. (Bottom) Intensities of the OH lines detected with NIRSpec. The 
uncertainties are only the fitting error. The impact of calibration effects is not taken into account and the associated uncertainties can be as high as 20%. 
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A (um) | Transition | Eup/kg (K) | 1 (x 1074 erg cm ? s7! sr!) 
17.085 | 0-0 S(1) 1015 1.55 + 0.19 
12.279 | 0-0 S(2) 1682 2.26 + 0.30 
9.665 0-0 S(3) 2504 12.85 + 0.68 
8.025 0-0 S(4) 3474 8.70 + 2.03 
6.910 0-0 S(5) 4586 20.27 £0.12 
6.109 0-0 S(6) 5830 5.30 + 0.12 
5.511 0-0 S(7) 7197 10.01 + 0.20 
5.053 0-0 S(8) 8677 1.96 + 0.15 
4.695 0-0 S(9) 10262 3.31 + 0.41 
4.410 | 0-0 S(10) 11940 0.52 + 0.33 
4.181 0-0 S(11) 13703 1.05 + 0.16 
3.846 | 0-0 S(13) 17444 0.47 + 0.07 
3.724 | 0-0S(14) 19403 0.24 + 0.13 
3.626 | 0-0 S(15) 21412 0.15 0.09 
2.627 1-0 O(2) 5987 9.89 + 3.15 
2.407 1-0 Q(1) 6149 35.70 + 0.62 
2.223 1-0 S(0) 6472 7.06 + 0.22 
2.122 1-0 S(1) 6952 25.03 + 0.36 
2.034 1-0 S(2) 7585 7.04 + 0.20 
1.958 1-0 S(3) 8365 14.23 + 0.42 
1.892 1-0 S(4) 9287 2.87 £0.13 
1.836 1-0 S(5) 10342 4.16+0.19 
1.788 1-0 S(6) 11522 0.42 + 0.17 
1.748 1-0 S(7) 12818 1.06 + 0.29 
2.786 2-1 O(2) 11636 0.65 + 0.13 
2.974 2-1 O(3) 11790 0.64 + 0.08 
2.356 2-1 S(O) 12095 0.20 + 0.08 
2.248 2-1 S(1) 12550 0.52 + 0.18 
2.073 2-1 S(3) 13891 0.47 + 0.14 
2.604 2-1 Q(5) 13891 0.27 + 0.15 
2.654 2-1 Q(7) 15763 0.19 + 0.08 
1.853 2-1 S(7) 18107 0.28 + 0.19 
1.233 3-1 S(1) 17819 0.72 + 0.10 
2.201 3-2 S(3) 19086 0.24 0.17 
2.066 3-2 S(5) 20857 0.16 + 0.02 
3.376 4-3 O(3) 22080 0.10 + 0.03 
1.242 4-2 S(4) 24734 0.15 + 0.09 


Extended Data Table 2. Intensities of the H3 lines detected at the bright spot with the OFF position subtracted. The uncertainties are only the fitting error. 
The impact of calibration effects is not taken into account and the associated uncertainties can be as high as 20%. 
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